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We have recently considered cosmologies in which the Universal scale factor varies as a power 
of the age of the Universe and concluded that they cannot satisfy the observational constraints 
on the present age, the magnitude- redshift relation for SN la, and the primordial element (D, 
3 He, 4 He, and 7 Li) abundances. This claim has been challenged in a proposal that suggested 
a high baryon density model (J7b^ 2 — 0.3) with an expansion factor varing linearly with time 
could be consistent with the observed abundance of primoridal helium-4, while satisfying the age 
and magnitude-redshift constraints. In this paper we further explore primordial nucleosynthesis in 
generic power-law cosmologies, including the linear case, concluding that models selected to satisfy 
the other observational constraints are incapable of accounting for all the light element abundances. 
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We have studied a class of cosmological models in which the Universal scale factor grows as a power of the age of the 
Universe (a oc t a ) and concluded that such models are not viable since constraints on the present age of the Universe 
and from the magnitude-redshift relation favor a — 1.0 ± 0.2, while those from the abundances of the light elements 
\ produced during primordial nucleosynthesis require that a lie in a very narrow range around 0.55 0. Successful 

■ primordial nucleosynthesis provides a very stringent constraint, requiring that a viable model simultaneously account 
^0 ' for the observationally inferred primordial abundances of deuterium, helium-3, helium-4 and lithium-7. For example, 

■ if the nucleosynthesis constraint is satisfied, the present Universe would be very young; to = 7.7 Gyr for a Hubble 
' parameter Ho = 70 kms -1 Mpc -1 (or, requiring Ho < 54 kms _1 Mpc _1 for to > 10 Gyr). 

Recently, Sethi et al. |2| noted that cosmologies where the scale factor grows linearly with time may produce the 
correct amount of 4 He provided that the Universal baryon fraction is sufficiently large. At first this result might 
seem counter-intuitive since such a Universe would have been very old at the time of Big Bang Nucleosynthesis 
(BBN) suggesting that all neutrons have decayed and are unavailable to be incorporated in 4 He. In fact, as Sethi et 
al. correctly pointed out, the expansion rate is so slow that the weak reactions remain in equilibrium sufficiently long 
to permit a "simmering" synthesis of the required amount of 4 He. However, such an old Universe also leaves more time 
to burn away D and 3 He so that no astrophysically significant amounts can survive. The observations of deuterium in 
high- redshift, low-metallicity QSO absorbers 0), the observations of lithium in very old, very metal-poor halo stars 
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(the "Spite plateau") ||, and those of helium in low-metallicity extragalactic Hn regions || require an internally 
consistent primordial origin. The Sethi et al. claim that deuterium could have a non-primordial origin is without basis 
as shown long ago by Epstein, Lattimer & Schramm ||. Nevertheless, the paper of Sethi et al. prompted us to 
reinvestigate primordial nucleosynthesis in those power-law cosmologies which may produce "interesting" amounts of 
He so as to study the predicted yields for D, 3 He, and 7 Li . 
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II. NUCLEOSYNTHESIS IN POWER-LAW COSMOLOGIES 

Preliminaries. For a power-law cosmology it is assumed that the scale factor varies as a power of the age independent 
of the cosmological epoch: 

a/a = (t/t ) a = (l + z)-\ (1) 

where the subscript '0' refers throughout to the present time and 'z' is the redshift. We may relate the present cosmic 
background radiation (CBR) temperature to that at any earlier epoch by T = (1 + z)/3Tq, where j3 < 1 accounts 
for any entropy production. For the models we consider, (3—1 after electron-positron annihilation. The Hubble 
parameter is then given by 
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The second equality should be read with the understanding that it is not valid during the epoch of electron-positron 
annhilation due to the non-adiabatic nature of annhilations. Power-law cosmologies with large a share the common 
feature that the slow Universal expansion rate permits neutrinos to remain in equilibrium until after electron-positron 
annihilation has ended so that neutrino and photon temperatures remain equal. In this case the entropy factor for 
T > m e /3 in eq. |^ is, 

13 = (29/43) 1/3 (3) 

in contrast to the standard Big Bang nucleosynthesis (SBBN) value of (4/11) 1 / 3 . As a increases, the expansion rate 
at a fixed temperature decreases due to the dominant effect of the 1/a power. Another useful way to view this 
is that at a fixed temperature, a power-law Universe with a larger a is older. As a consequence of the decreasing 
expansion rate, the reactions remain in equilibrium longer. In particular, as pointed out by Sethi et al. [0 for the 
linear expansion model (a = 1), the weak interactions remain in equilibrium to much lower temperatures than in the 
SBBN scenario, allowing neutrons and protons to maintain equilibrium at temperatures below 100 keV, as can be 
seen in Fig. |l|. As is evident from Fig. |l|, the 4 He production rate below about 0.4 MeV is too slow to maintain 
nuclear statistical equilibrium. However, the presence of neutrons in equilibrium and the enormous amount of time 
available for nucleosynthesis during neutron-proton equilibrium (compared to SBBN) make it possible to build up a 
significant abundance of 4 He 0. 
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FIG. 1. Comparison of nucleosynthesis in the linear expansion model (heavy curves) for the case of Yp = 0.24 with the 
predictions of nuclear statistical equilibrium (lighter curves). The solid curves are for the 4 He mass fraction Yp, while the 
dashed curves show the evolution of the ratio of neutrons to protons. 

The above discussion is not restricted to a = 1, but applies for all values of a which are sufficiently large (so that 
the expansion rate is sufficiently small) to allow neutrons to stay in equilibrium long enough to enable synthesis of 
4 He in sufficient amounts, as we show in Fig. |^. Although we explore a larger range in a in this paper, we present 
detailed results for 0.75 < a < 1.25, a range consistent with the age and expansion rate of the Universe, and we check 
these results for consistency with independent (i.e., non-BBN) constraints on the baryon density. The iso-abundance 
contours in Fig. 2 show clearly that as a decreases towards 0.75, a larger baryon density is required to produce the 
same abundance of 4 He. For example, although Yp = 0.24 can be synthesized in the a = 0.75 model, the density of 
baryons required is very large: Slp,^ 2 — 20. These 'large f2p,' models are constrained by dynamical estimates of the 
mass density, an issue we discuss later. 
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FIG. 2. Iso-abundance contours of the 4 He mass fraction (Yp) in the baryon density - a plane. The shaded band corresponds 
to helium abundances in the range 0.22 < Yp < 0.26. 

Helium-^ abundance. In an earlier study Q, we showed that there is a very small region, centered on a = 0.55, for 
which the light elements can be produced in abundances similar to those predicted by SBBN. But, this small window 
is closed by the SNIa magnitude-redshift data |^|. Here we are concerned with larger values of a and, correspondingly, 
larger baryon-to-photon ratios (r/). First we consider the nucleosynthesis of 4 He in these models. Figure^ shows the 
connection between the baryon density (fish 2 — fyio/273, where ijio = 10 10 npf /n 7 ) and a set by the requirement that 
the primordial helium mass fraction lie in the generous range 0.22 < Yp < 0.26. We have included in Fig. || the region 
investigated in Ref. a < 0.6, as well. To understand the features in Fig. g, we need to isolate the important factors 
controlling the synthesis of helium. In SBBN, the 4 He abundance is essentially controlled by the number density ratio 
of neutrons to protons (n/p) at the start of nucleosynthesis (T = Tbbn ~ 80 keV). This ratio in turn is determined 
by (1) the n—p ratio at 'freeze-out' (T — Tf) of the neutron-proton interconversion rates which may be approximated 
by {n/p)f = cxp(—Q/Tf), where Q = 1.293 McV is the neutron-proton mass difference and, (2) the time available 
for neutrons to decay after freeze-out, Atd = t(TsBN) — t(Tf). In contrast, for power-law cosmologies another factor 
comes into play - the time available for nucleosynthesis, AtsBN, before the nuclear reactions freeze-out. For larger 
a, the expansion rate of the Universe (at fixed temperature) is smaller and the Universe is older. Hence, for larger 
a neutrons remain in equilibrium longer and the freeze-out temperature (Tf) is smaller, so that (n/p)t is smaller. 
However, the effect of the increase in AtsBN as a increases dominates that due to the change in Tf. For a — 0.50 the 
freeze-out temperature is around 4 MeV whereas for a — 0.55, Tf ~ 1 MeV which implies a decrease in (n/p)f by a 
factor of about 2.5. On the other hand, the age of the Universe at T = 10 keV (about the temperature when SBBN 
ends) is a factor of 25 larger for a = 0.55 relative to that for a = 0.50. Thus, for the same rj, increasing a from 0.50 
to 0.55 has the effect of increasing the 4 He abundance because more time is available for nucleosynthesis. But, since 
decreasing the baryon density decreases the nuclear reaction rates leading to a decrease in 4 He, we may understand 
the trend of the smaller baryon density requirement as a increases from 0.50 to about 0.55, even though the decrease 
in Tf opposes this effect. The time-delay between 'freeze-out' and BBN, Atd, which has, until now, been much smaller 
than t„, becomes comparable to it at a ~ 0.55. Since a larger a results in an older Universe at a fixed temperature, 
Atd increases with a. Thus for a > 0.55, Yp is increasingly suppressed (exponentially) as a is increased. The only 
way to compensate for this is by increasing Tbbn (since Atd oc (Tbbn) 1 ^ ), which may be achieved by increasing 
the baryon density. But since Tbbn depends only logarithmically on the baryon density fiofl , this accounts for the 
exponential rise in the required value of fish 2 as a increases. This trend cannot continue indefinitely; the curve must 
turn over for reasons we describe below. 
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From Fig. [|, it is apparent that in the "large a" range, the required value of VL^h 2 decreases with increasing a. In 
our previous analysis Jlj] of 4 He nucleosynthesis which concentrated on a in the vicinity of 0.55, we implicitly assumed 
that the age of the Universe at T — Tf was not large enough for appreciable amounts of 4 He to have been built 
up. This assumption breaks down for large values of a and r\. Since D, 3 He and 3 H are not present in appreciable 
quantities, a large value of r\ is needed to boost the 4 He production rate. Now, the larger the value of a, the longer 
neutrons remain in equilibrium, thus allowing more 4 He to be slowly built up, with the neutrons incorporated in 4 He 
being replaced via p — > n reactions. Roughly speaking, the required value of 77 for a given a is set by the condition: 
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The effects of a on t(Tf) and n on dYp/di complement each other, giving rise to the trend shown by the 4 He 
iso-abundance curves in Fig. @ for a > 0.75. 
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FIG. 3. Evolution of the light element abundances as a function of the photon temperature in an a = 1 Universe. 

Light element abundances in the linear expansion model. We now turn to the production of deuterium and 3 Hc. 
For large a {e.g., a = 1), we expect the deuterium abundance to be insignificant since D can be efficiently burned to 
3 He during the long time available for nucleosynthesis. The mean lifetime of deuterium against destructive collisions 
with protons at a low temperature of 10 keV is around 3 days; at this temperature the a = 1 Universe is already 300 
years old! The fact that the timescales are so different allows us to derive analytical expressions for the deuterium, 
helium-3 and lithium-7 (beryllium-7) mass fractions (to be denoted by Xjj, X3 and X7 respectively). The generic 
equation for the rate of evolution of the mass fraction of nuclide "a" can be parameterized as 

— J^- = ^prod(a) — -Rdost(a)^a, (5) 

where "prod" and "dest" refer to the production and destruction rates of nuclide "a" . Given that the universe remains 
at the same temperature for a very long time (compared to the reaction time scales), it is not surprisingthat X a 
achieves its steady-state value at each temperature (for a detailed discussion in the context of SBBN, see jyj), 

„ _ -Rprod(a) 

a ~ ~R 7n\- W 
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We can write this explicitly for the simplest case - deuterium: 



X 



D 



pD 



(7) 



where the various Ts represent the relevant deuterium creation (n + p — > D + 7 and p + p — > D + e + + v) rates 
per target proton, and destruction (D(p, 7) 3 He and D(pj,p)n) rates per target deuterium. All of these rates can be 
obtained from Ref. jl5). Once the reaction rates become smaller than the universal expansion rate (say at some 
temperature T*), the abundances freeze out with values close to X a at the corresponding T*. This is illustrated in 
Fig. ^ which clearly shows that the steady-state solution works very well. We note here that the steady state (dotted) 
curves in Fig. [| are not independent analytic derivations, but use the abundances of the various nuclei as calculated 
by the numerical code. The figure intends to emphasize that nucleosynthesis in this (linear expansion) model can be 
well represented by the steady-state solutions in eq. ^. 
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FIG. 4. Comparison of the light element abundances with their steady state values as a function of the photon temperature 
in an a = 1 Universe. The dotted curves correspond to the equilibrium solution. 

In the expression for Xu (see eq. [j]), the n + p reaction term dominates until about 20 keV after which the p + p 
reaction makes the dominant contribution. The final deuterium abundance is thus determined by the weak pp reaction 
{p + p — > D + e+ + u), the effect of which can be seen in Fig. || as the very slow rise in X D between temperatures 
of 10 keV and 1 keV (at which point the D abundance freezes out). Since both 3 He and 7 Li freeze out much earlier, 
they do not get any significant boost from the weak pp reaction. 

From eq. Xu, and thus X3 ( 3 He is formed from D), are proportional to X n , the neutron abundance. One striking 
feature in Fig. || is the boost to the neutron abundance (and hence the abundances of D and 3 He) at temperatures 
around 40 keV. The effect is subtle and may be missed in BBN codes with a limited nuclear reaction network. The 
slow rate of expansion of the universe during nucleosynthesis facilitates the production of a relatively large "metal" 
(A> 8) abundance (X meta is — 3 x 10~ 7 ). In particular, 13 C is produced in these models through the chain: 12 C+p — ► 
13 N+7 and the subsequent beta-decay of 13 N. In this environment 13 C+ 4 He — > 16 0+n leads to the production of free 
neutrons. 

The mass-7 abundance is entirely due to the production of 7 Be through the reaction 3 Hc+ 4 Hc — > 7 Bc+7. 7 Bc decays 
to 7 Li by electron capture once the universe has cooled sufficiently to permit the formation of atoms. Once formed, 
it is difficult to destroy 7 Be at temperatures < 100 KeV. In contrast, 7 Li is very easily destroyed, specifically through 
its interaction with protons. Since the 7 Be production (and thus 7 Li) follows the evolution of the 3 He abundance, 
and there is very little destruction of 7 Be, 7 Li also benefits from the boost to the neutron abundance described in 
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the last paragraph. This has the effect of boosting the 7 Li abundance from 10 11 (if this source of neutrons were not 
included) to 1CT 9 . This is significant in that, at the level of a few parts in 10 10 (e.g., the primordial lithium 

abundance lies between these two estimates. 
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FIG. 5. Abundances of 3 He, D and 7 Li in power-law cosmologies for different values of the expansion index (a). The shaded 
bands correspond to helium abundances in the range 0.22 < Yp < 0.26. 

Light element abundances vs. a. Having explored BBN in the linear model (a = 1) it is now important to ask how 
these results depend on a. It is clear from Fig. ^ that nothing dramatically different happens as a changes; this is 
simply because the key physics remains the same. In preparing Fig. ^| we adjust the value of rj (baryon density) for 
each choice of a so that the primordial 4 He mass fraction lies between 22% and 26%. As a increases, the nuclei freeze 
out at lower temperatures since the expansion rate at the same temperature is lower for a larger a. The effect of this 
can be gauged by the behavior of Xd, X$ and X7 with respect to temperature, as given by eq. ^]. For deuterium this 
implies a small increase with a due to the pp (weak) reaction, which is also reflected in the behavior of 3 He for a > 1. 
The fall of 3 He with increasing a for a < 1 is due to larger destruction of 3 He because of the increase in the time 
available for the nuclear reactions. As already mentioned, the abundance of 7 Be depends critically on the evolution of 
the 3 He abundance; so while the mass-7 ( 7 Be) abundance increases appreciably with the increase in baryon density, 
it is relatively unaffected by a change in a. 

Note that in those power-law models which can simultaneously reproduce an acceptable 4 He abundance along 
with a consistent age and expansion rate, the corresponding baryon density must be very large, 0.04 < ^l^ti 2 < 6.4 
(11 < j]i < 1750; see Fig. 2). Most - if not all - of this range is far too large for consistency with independent (non- 
BBN) estimates of the universal density of baryons (rj < 7.4 Jl6|) or, for that matter, the total matter density ||. 
Conservatively, clusters limit the total (gravitating) matter density to Qm 0.4 so that if there were no non-baryonic 
dark matter, 0,-Qh 2 < 0.2 (77 < 54) for h ~ 0.7. However, if the X-ray emission from clusters is used to estimate the 
cluster baryon fraction (see fll7|]), the universal baryon density should be smaller than this very conservative estimate 
by a factor of 7-8 (consistent with the upper bound from the baryon inventory of Fukugita, Hogan and Peebles jTJJ ) . 
Thus power-law cosmologies constrained to reproduce He (only), an acceptable age and magnitude-redshift relation, 
and an acceptable baryon density, must have a restricted to a very narrow range: 1 < a < 1.2. Furthermore, the 
baryon density in even this restricted range is large when compared with estimates [ p7[ of the baryon density from 
cluster X-rays. Finally, for a in the narrow range of 1 < a < 1.2 and 0.22 < Yp < 0.26, the other light element 
abundances are restricted to 7 Li/H > 10~ 9 , 3 He/H < 3 x 10~ 13 , and D/H < 3 x 10~ 18 . For deuterium and helium-3 
this is in very strong disagreement (by 8 to 13 orders of magnitude!) with observational data (for a review see |jl8||). 
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Although the predicted 7 Li abundance is comparable to that observed in the solar system, the local ISM, and in 
Pop i stars, it is larger than the primordial abundance inferred from the Popn halo stars f3lJl8f, and marginally 
inconsistent with the observations of lithium in the ISM of the LMC . 



III. CONCLUSIONS 



In response to the claim J2J that a power-law Universe expanding linearly with time could be consistent with the 
constraints on BBN, we have reexamined these models. Although it is true that observationally consistent amounts 
of 4 He can be produced in these models, this is not the case for the other light elements D, 3 He, 7 Li. Furthermore, 
consistency with 4 He at a = 1 requires a very high baryon density (75 < ijiq < 86 or 0.27 < Vl-^h 2 < 0.32), 
inconsistent with non-BBN estimates of the universal baryon density and, even with the total mass density. We 
have also investigated BBN in power-law cosmologies with a > 1 and have confirmed that although the correct 
4 He abundance can be produced, the yields of the other light elements D, 3 He, and 7 Li are inconsistent with their 
inferred primordial abundances. In general, power-law cosmologies are unable to account simultaneously for the early 
evolution of the Universe (BBN) (which requires a ~ 0.55) and for its presently observed expansion (which requires 
a = l±0.2). 
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